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Three-Body Collective Excitations in the Superconducting Phase of 

MgB2 

Naoum KarchevQ 
Department of Physics, University of Sofia, 1164 Sofia, Bulgaria 

It is shown that coherent behavior of Cooper's pairs and electrons in the two-band superconductors 
MgB2 is possible as a result of a strong two-phonon-electron coupling. The spin-i triples with 
zero angular momentum are made up of three spin-i fermions with charge e. They are gapped 
Fermi excitations with gap induced by the gaps of the single fermions. Their contribution to the 
thermodynamics of the MgB2 superconductors is considered. 
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I INTRODUCTION 

The three-body bound states are everywhere in the 
physics, but the study of their role in collective condensed 
matter behavior is still limited. Phase transitions driven 
by an instability in a three-fermion channel have been 
explored. The anomalous three-body scattering ampli- 
tude is used as an ansatz to develop a mechanism of 
odd- frequency superconductivity 0. Unhke a BCS the- 
ory, the authors involve a cooperative pairing of electrons 
and spins. The three-body bound state formation leads 
to gap which is an odd function of frequency. The appli- 
cation to heavy fermion superconductors is discussed. 

In the present paper the interest in this topic is in- 
spired by the anomalous superconducting properties of 
MgB2 |3| • The superconductivity in magnesium diboride 
is s wave, mediated by electron-phonon coupling. It dif- 
fers from ordinary metallic superconductors in several 
ways. Scanning tunneling microscopy Q and point con- 
tact studies revealed double-peaked spectra at low 
temperature that were interpreted as evidence for two 
gap superconductivity. Ab initio calculations suggest 
that multiple gaps are a consequence of the coupling 
of distinct electronic bands |^ |g . MgB2 has strongly 
anisotropic Fermi surface of four separate sheets that 
are grouped into two-dimensional a bands and three- 
dimensional TT bands. The different energy gaps are as- 
sociated with the TT and a sheets. The key quantity 
in phonon-mediated superconductivity is the Eliashberg 
function a^(a;)_F'(w), where F{uj) is phonon density of 
states and a{uj) is the electron-phonon coupling averaged 
over directions in k space. The electron tunneling spec- 
troscopy is used to determine the three distinct Eliash- 
berg functions (a^F)7r-7r, {a'^F)cr-cr, and (a^F)^^-^- 

The double-gap structure is used to explain some of 
the unusual physical properties of MgB2, such as the 
rapid rise of the specific heat coefficient C/T j^, tun- 
neling 5] and upper critical field anisotropy Quite 
different methods of theoretical investigation, includ- 
ing weak-coupling two-band BCS theory ^3|> Eliash- 
berg strong-coupling formalism 6], and strong-coupling 
density-functional technique with explicit account for the 
Coulomb repulsion lead to astonishingly identical 
curves for specific heat, as a function of temperature. 



The calculations reproduce the different slopes, above 
O.bTc and below 0.25rc, which apparently result from 
the existence of two gaps, but can not explain the shoul- 
der between them ■ Unexpected are the effects of 
Mg substitution by Al on the specific heat. The changes 
are in rather poor agreement with those predicted by tak- 
ing into account changes in the electronic and photonic 
structure only [l3l |. 

A large B isotope effect is another argument in favor 
of phonon-mediated pairing 0, ^] . The isotope coeffi- 
cient a is defined by the relation Tc oc M~°', where M is 
the mass of the element. In BCS theory a = 0.5, and for 
metals like Hg, Pb and Zn the coefficient is found exper- 
imentally to be close to 0.5. The isotope coefficient for 
MgB2 is a « 0.32. The density-functional calculations 
of the phonon spectrum and electron-phonon coupling in 
MgB2 predict that in this compound, phonon modes of 
Boron oscillations may have relatively high frequencies, 
and that nonlinear coupling via two-phonon exchange is 
comparable to or even larger than the linear coupling 
Both effects may contribute to the anomalous iso- 
tope effect coefficient [l^, and to the significant increas- 
ing of the critical temperature Tc — 39K. It is thought 
that Tc of MgB2 probably represents the upper limit of 
the phonon mediated superconductivity. 

Motivated by the theoretical and experimental findings 
I consider theory of two-band superconductors [TtI 
with two-phonon electron interaction. The main goal is 
to study the formation of coherent behavior of Cooper's 
pairs and electrons in the two-band superconductors 
MgB2, as a result of a strong two-phonon-electron cou- 
pling. It is shown that spin-i triples with zero angu- 
lar momentum, made up of three spin-i fermions with 
charge e, are possible. They are gapped Fermi excitations 
with gap induced by the gaps of the single fermions. Ef- 
fectively one can represent them as gapped fermions and 
to write an effective action. This enable to calculate the 
contribution of the triples to the thermodynamics of the 
MgB2 superconductors. To reproduce the shoulder in 
the specific heat as a function of temperature, one has to 
choose the gaps of the triples larger than the gaps of the 
incipient electrons. 

The paper is organized as follows. In Sec. II the for- 
mation of triples and their contribution to the specific 



2 



heat is explored. A summary in Sec. IV concludes this 
paper. The tunneling experiments are debated as a most 
direct way to observe the triples experimentally. 



II TRIPLES IN TWO-BAND SUPERCON- 
DUCTORS WITH NONLINEAR ELECTRON- 
PHONON COUPLING 

I consider theory of two-band superconductors 0, 
with two-phonon electron interaction. An important con- 
sequence of this interaction is the effective six-fermion 
interaction. One can obtain it from the triangle diagram 
(Fi^l with three phonon lines (undulating lines). There 



where /3 is the inverse temperature and e^(V) is the dis- 
persion of band i fermion. The Hamiltonian iJint('7") is a 
sum of the BCS Hamiltonian jSJl and six-fcrmion Hamil- 
tonian J^l- It is obtained from Eqs.(^) and ||2J) replacing 
the operators with Grassmann functions. 

Let us introduce two spin-i fermi collective fields 
(triples) Cfcr(T, x) (Ci(T, x)) by means of the Hubbard- 
Stratanovich transformation of the six-fermion term 

e-"f^ = j i?A^(C+C)exp{- j d^xJ2^i [CtMQ<rix)+ 



(6) 



CfT ix)cji (x)ctg„ {x)Cla (x) + {x)C-la ix)cil {x)cl^^ ix) } 




FIG. 1: Triangle diagram with three phonon lines (undulating 
lines) which results from the two-phonon exchange 

are two fermion species in the theory, therefore the low 
frequency and low momenta limit of the diagram leads to 
an effective local six-fermion term. The effective Hamil- 
tonian of the six-fermion interaction has the form 



(1) 



- / d^xc[^ (x)c+ (x)c+^^ (x)c_£<, (x)q4 (x)q| (x) 
No- 
where c^(x) and CfCT(x) are creation and annihilation op- 
erators for fermions, with orbital index £ (£ = 1, 2, —I = 
2, 1) and spin projection a (a =t, J.) [l9l |. 

The BCS Hamiltonian of the theory of two-band su- 
perconductivity is 0, 

Hbcs = -J2 9e j d^xcj^{x)cj^{x.)cii{x)c^{x) (2) 

-33 Jd^xJ2 4-\ (x)c-a Wc^^T (x) 

The partition function can be written as a path inte- 
gral over the Grassmann functions of the Matsubara time 
T c^(r,x) and Qcr(T,x) 



(3) 



Z{(3) = j Dfi{c+c)e- 
The action is given by the expressions 



S ^ So + f dTHi^tir) 



(4) 



/3 



So^JdrJ d3xEct(^,x)[a, + e,(V)]c,,(r,x) (5) 
i 



where x = (t, x) and J dr J cPx — J d'^x. 



The effective action of the triples is defined by the 
equality 



e-SeijCC+.O^expj-l d'xJ^X,Cl{^)Q4x)^ (7) 

exp{~ J d*x^Xe ^cj^{x)cj^{x)c'ti^{x)Cea{x) + 
0^(2;)c-f^(a;)Qx(a;)Q|(x)]} > 



< 



with 



<Q>= J Dfi{c+,c)Qe-^'>- 



(8) 



The quadratic part of the effective action Seff (C^ , C) has 
the form 



S. 



eff 



d'^xd'^y 



<:t{x)Ilii,{x~y)Ce'a'{y)+{9) 



Cl{x)Ki'ix ~ y)C^Ay) + CA^)^.Ax - y)Ci'a'{y) . 

where 

Hf^,(a; ~y)^ 5ij'5aa'\i5'^{x - y) - 
XiXt < c-ia{x)cii{x)ci'\{x)cj,^{y)cj, ^{y)cti,^,{y) >(10) 

KUx-y) = 



2 

X^Xii 



< c-i^{x)cii{x)c^{x)c-£>cr'iy)ci'iiy)ci'i;{y) >(11) 

a' 

^aa'{x-y) = 

< c+{x)4^{x)cte^ix)cj,^iy)cj,^iy)c'te,„,iy) :=(12) 



In theory with Hamiltonian Hbcs © the off-diagonal 

elements are zero: H"^, = 0, E"^, = and S^^./ = if 
£ ^ £' . The diagonal functions are calculated in leading 
order represented by the diagrams (Fig2). In the case of 
H^^^, (Fig2a), two lines of the diagram (solid lines) cor- 
respond to the normal Green function of fermions with 
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one and just the same band-index, while the third hne 
(dashed Hne) corresponds to the normal Green function 
of fermion with different band-index. The diagrams for 

E^^^, and S^^, (Fig2b) have two hues (solid hues) corre- 
sponding to the anomalous Green functions of fermions 
with equal band-index and one line (dashed line) which 
corresponds to the anomalous Green function of fermion 
with different band-index. 



(27r)353(p^+P2+P3-p) 



n (27r)3 4£;_,p^^,p^£;,p^(i;_,p^ + Sfp^ + S,pJ 
£-fp^Ei-p^€i-p^ + e-fPjQp^ii^fPg] 



(18) 



a) 




££(P) = 1- 



(27r)353(pi+p2+P3-p) 



b) 




^ (27r)3 4£'_£p^i?£p^£'£p^(i?_fp^ + Sfp^ + -B^Pg) 
[e-fp^e^p.e^Pg + £-tp,Ei^^Et^^ + 
E-ip^Sip^Eep^ + E^ep^Eip^Eep^] 



(19) 



FIG. 2: Two loop diagrams which represent leading order of: 
a) n«„b) E«, and sit, 

The system of gaps' equations in the theory of two- 
band superconductivity has the form 0, 0| 



Mp) = 

d?p^ A^A-^(27r)3,53(p^-Hp2-Hp3-p) 
2 J (27r)3 ^E^tp,Etp^Eip,^{E^ip^ + Eep^ + Eep^ 



(20) 



Ai = 



91 



tanh(§i;ip 
Eip 



■Ai 



.93 



(Tp 



tanh(fi;2p) 



E' 



A, 



The equations 



2p 



(13) 



£i{Pfi) = 



(21) 



53 



3 tanh(f£;ip) g2 f 3 tanh(f£;2p define the Fermi surface of the triples. Using the approx- 
A2 = — I d p ^^7— Ai + I d. p ^^7"^^; ^2 imate expressions for Zi{p) and Sf(p) 



Ell 



E2p) 



with Egp = \J ^Ip + l^f |j where etp is ^-band fermion 

dispersion. The system of gaps' equations shows that 
the gaps Ai can be chosen real. Then the normal and 
anomalous Green functions in two-band theory, calcu- 
lated in the standard way, have the form [23 



5^^(w,p) = Sl^{uj,p) = 
F{uj,p)^ F^{lo,p) 



£ip 



^2 + + A2 



I calculate the diagrams in low-frequency limit 
n« ,(cc>,p) = 5„„, [-iujZ-\p) + e,(p)] 



S«(0,p) = S"(0,p) = E,(p 



(14) 
(15) 

(16) 
(17) 



Zt (p) ~ Z, [pf, )^Ze, Mp)-^i (Pfi ) , 
and re-scaling the triples' fields 

z'Hl{^,p)^ct{^,p), 

Z~iCeaiuJ,p) -> 0,7 (w,p) 
one obtains the effective action of the triples 

P 

Seff = [drY, {(tir, P) {Or + iiip)) Qa{T, p) + 
P 

Qi [c^ (t, p)C^ (t, -p) + Cn (t, -p)Cn p) } 



(22) 

(23) 
(24) 

(25) 



The result is 



with £f(p) — Z^ieip) ^^"^ — ZgYii{pjf). Near the 

Fermi surface £i{p) — ^{p — pfi), where rhi are triples' 
masses. 

The effective action shows that triples are spin-| 
fermi excitations with gap gi induced by the gaps of the 
single fermion excitations (|20|l . Next one diagonalizes the 
effective Hamiltonian using a Bogoliubov transformation, 
and rewrites it in terms of Bogoliubov excitations with 
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9l 



This enable to calcu- 



dispersion Eeip) = \/ ^'j{p) 
late the contribution of triples to the thermodynamics 
of superconductors. In particular the low temperature 
behavior of the heat capacity is 



E 



mipft 



2nA° 
~T3~ 



2TTPf _Pl_ 

— —e T 

T3 



(26) 

where the first terms come from the single fcrmion con- 
tribution, while the other terms are triples' contribu- 
tion. Recent experiments, including photoemission 
and tunnehng experiments S H 

, show that the ratio 
of the single fermions' gaps is 2.6 < A2/A1 < 3.5. The 
low temperature behavior of the heat capacity coefficient 
— — 7 is depicted as a function of in Fig. 3 

= 4, ^ = 5 and mipji ~ m2P/2 ^ 
The curve is in a good agreement with 



mip fi\'27r 

for 1^ = 



Ai 



the experimental one for MgB2 
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FIG. 3: The heat capacity coefficient — 7=77? 

function of ^ for ^ = 3, = 4, = 5 and m\Pfi 
m2P/2 — rhipfi ~ rri2Pf2- 



III SUMMARY 

In summary, new type of excitations, triples- made 
up of three spin-i fermions, are predicted in theory of 
two-band superconductivity with non-linear two-phonon- 
fermion interaction. They can be thought of as a bound- 
state of s-type Cooper pair of fermions from one of the 
bands and fermion from the other one with zero angu- 
lar momentum. The triples are gapped excitations with 
gap induced by the single fermion gaps. It is important 
to underline that the triples result from the nonlinear 
two-phonon-electron interaction. One can consider exci- 



tations made up of more than three spin-^ fermions while 
studying diagrams with more than three vertexes. Due to 
the Pauli principle these excitations have non-zero angu- 
lar momentum. The specific form of the fermion interac- 
tions determine the symmetry of these excitations. They 
in turn result from the non-linear two-phonon-fermion 
interaction. Therefore the knowledge of the non-linear 
phonon interaction is crucial for the development of a the- 
ory of excitations made up of more than three fermions. 
The same is true for triples made up of fermions from one 
band. They have non-zero angular momentum too, and 
the symmetry of the triples results from the non-linear 
two-phonon-fermion interaction. 

The six-fermion interaction ^ can be alternatively ob- 
tained from the linear phonon-electron coupling. The 
leading order diagram is depicted in Fig4. Since my in- 




FIG. 4: Six-fermion interaction obtained from the linear 
phonon-electron coupling 

vestigation is based on the assumption that nonlinear 
coupling via two-phonon exchange is comparable to or 
even larger than the linear coupling ^sl, 'l(^ , the contri- 
bution of the linear phonon-electron coupling (Fig4) is 
small perturbation to the contribution of the nonlinear 
phonon-electron interaction and one can drop it. 

It is easy to see from the formula for the heat capacity 
(|25|l that the contribution of the triples is decisive to 
explain the shoulder like part of the curve. 

The most promising way to observe the triples in 
MgB2 is by tunneling experiments. It is evident that 
the contribution of the triples to the tunneling current is 
with much smaller weight than those of the single elec- 
trons. Therefore, we can observe the triples only at very 
low temperature. If one considers the differential con- 
ductance of tunnel junction with counter-electrodes in 
superconducting state, for example In, or Ph as in 0, 
two new peaks should emerge when the temperature de- 
creases. Tunneling experiments achieved bellow \K can 
answer the question about the existence or nonexistence 
of the triples. 
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